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We perform an amplitude analysis of the decays B° — > 07^1(1430)°, if* (892)°, and </>(if7r)g_ wavc 
with a sample of about 384 million BB pairs recorded with the BABAR detector. The fractions of 
longitudinal polarization /l of the vector-tensor and vector-vector decay modes are measured to be 
0.8531 

o!o69 ±0.036 and 0.506 ±0.040 ±0.015, respectively. Overall, twelve parameters are measured 
for the vector- vector decay and seven parameters for the vector-tensor decay, including the branching 
fractions and parameters sensitive to CP- violation. 

PACS numbers: 13.25.Hw, 13.88.+e, 11.30.Er 



The interest in the polarization and CP-asymmetry 
measurements in B — > <j)K* decays is motivated by 
their potential sensitivity to physics beyond the standard 
model in the b — > s transition, shown in Fig. [1] (a) [![. 
The polarization measurements of B meson decays re- 
veal both strong and weak interaction dynamics and are 
discussed in a recent review 0, The large fraction 
of transverse polarization in the B — > 0if*(892) decay 
measured by BABAR |4| and by Belle [5] indicates a sig- 
nificant departure from the naive expectation of predom- 
inant longitudinal polarization. This suggests other con- 
tributions to the decay amplitude, previously neglected, 
either within or beyond the standard model [6| . 

We now extend our investigation of the polarization 
puzzle with an amplitude analysis of the vector-tensor 
B° — > 4>K2 (1430)° decay. We also measure vector- vector 
B° -> (j)K* (892)° and vector-scalar B° -> (/)(Kir)f decay 
amplitudes, where (Ktt)q° is the J p = + Ktt compo- 
nent. We use the dependence on the Ktt invariant mass 
of the interference between the J p = + and l - or 2 + 
components 0, Q to resolve the discrete ambiguity in the 
determination of the strong and weak phases otherwise 
present in the B° -> c/>K* (892)° analysis and to 

provide new measurements of the strong and weak phases 
relative to the vector-scalar decay amplitude. 

The angular distribution of the B — > cf>K* decay can 
be expressed as a function of Hi — cos 0i and $ shown 
in Fig. [1] (b) . Here i is the angle between the direction 
of the K meson from the K* — > Ktt [6{) or cj> — > KK 
(82) and the direction opposite the B in the K* or <j> rest 
frame, and <E> is the angle between the decay planes of 
the two systems. The differential decay width has seven 
complex amplitudes Aj\ corresponding to the spin of the 
Ktt system J and the helicity A — or ±1: 



We analyze 'i?' 



d 3 T 



dHidH 2 d<I> 



oc 



Y,Aj X Y}{H u *)Y{*(-H2,0) 



(1) 



where Yj are the spherical harmonics with J = 2 for 
#|(1430), J = 1 for K*(892), and J = for {Ktt)* . 
We can reparameterize the amplitudes with the index J 
suppressed as A® and A±i — (A\\ ± A±)/y2. 



K*° -> (K+K-^K^) candi- 
dates using data collected with the BABAR detector [9j at 
the PEP-II e+e" collider. A sample of 383.6±4.2 million 
T(4S) — > BB events was recorded at the center-of-mass 
energy s/s — 10.58 GeV. Charged-particle momenta are 
measured in a tracking system consisting of a silicon ver- 
tex tracker with five double-sided layers and a 40-layer 
drift chamber, both within the 1.5-T magnetic field of a 
solenoid. Charged-particle identification is provided by 
measurements of the energy loss in the tracking devices 
and by a ring-imaging Cherenkov detector. 

We use two kinematic variables: AE = {EiEb — Pi ■ 
p s -s/2)/Vs and m ES = [{a/2 ± p 4 • Pb) 2 /E? - pl] 1/2 , 
where (Ei,pi) is the e + e~ beam four-momentum, and 
(Eb, Pb) is the four-momentum of the B candidate. We 
require \AE\ < 0.1 GeV and m E s > 5.25 GeV. The re- 
quirements on the invariant masses are 0.99 < < 
1.05 GeV and 0.75 < mx-ix < 1-05 GeV (lower ma-* range) 
or 1.13 < ttikti < 1.53 GeV (higher range). 

To reject the dominant e + e~ — > quark-antiquark back- 
ground, we use variables calculated in the center-of-mass 
frame. We require |cos6>t| < 0.8, where 9t is the an- 
gle between the i?-candidate thrust axis and that of the 
rest of the event. We construct a Fisher discriminant, T, 
that combines the polar angles of the B-momentum vec- 
tor and the B-candidate thrust axis with respect to the 
beam axis, and the two Legendre moments Lo and L2 of 
the energy flow around the B-candidate thrust axis • 
We remove signal candidates that have decay products 




0i A K* 



d 



(a) 



d 





FIG. 1: (a) Feynman diagram describing the B° — > 4>K*° 
decay; (b) definition of decay angles given in the rest frames 
of the decaying parents. 
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with invariant mass within 12 MeV of the nominal mass 



values for D^r or 



far . In about 5% of events, 



more than one candidate is reconstructed, and we select 
the one whose four-track vertex has the lowest x 2 . We 
define the flavor sign Q to be the charge of the pion. 

We use an unbinned, extended maximum-likelihood 
fit Q to extract the event yields rij, flavor asymmetries 
Aj , and the probability density function (PDF) parame- 
ters, denoted by C for the polarization parameters and £ 
for the remaining parameters. The data model has five 
event categories j: B — > <f>(Kn)j=o,i,2, B —> f Q (980)K*, 
and combinatorial background. The combinatorial back- 
ground PDF is found to account well for both the domi- 
nant quark-antiquark background and the random tracks 
from the B decays. The likelihood £j for each candidate 
i is defined as = V • k n 1 - Vjfai; £; £), where each of 
the Vj is the PDF for variables = {Hi, Hi, TOKin 
m K x, AE, rriEs, J~, Q}- The flavor index k corresponds 
to the value of Q, that is Vj = Vj x S^q- 

We define rij = 71+ + nj and Aj = (rij — rij)/ (rij + 
rij). The polarization parameters, with the index J 
suppressed, are defined as Jl = \A n \ 2 /Y,\A\\ 2 , f± = 
|Al|7£|A a | 2 , 0|| = arg(A||/A)), and <f> x = &xg(A±/A ). 
We allow for CP-violating differences between the B° 
(Q = +1) and B° (Q = —1) decay amplitudes (A and 
A) and incorporate them via the replacements /l — > 
f L x(l+A° CP xQ), / ± -» f±x(l+A^ P xQ),<j> {] -» (0H- 
A0|| x Q), and ± -» (</>j_ + tt/2 + (A^ + tt/2) x Q) @j. 

The PDF Vj(y^i, C; £) for a given candidate i is a joint 
PDF for the helicity angles, resonance masses, and Q, 
and the product of the PDFs for each of the remaining 
variables. The helicity part of the exclusive B decay PDF 
is the ideal angular distribution from Eq. JT]), where the 
amplitudes Aj\ are expressed in terms of the polarization 
parameters multiplied by an empirically-determined 
acceptance function Q(Hi,H2,$) = Gi(Hi) x Qi^Hi)- 
A relativistic J-spin Breit-Wigner amplitude parameter- 
ization is used for the resonance mass 0, [ll| , except for 
the (Ktt)^P niK-ir amplitude parameterized with the LASS 



function [7[. The latter includes the Kq (1430)° resonance 
together with a nonresonant component. 

The interference between the J = 1 or 2 and the S- 
wave (Kir) contributions is modeled with the three terms 
21Ze(Aj\Ag ) in Eq. with the four-dimensional an- 
gular and rriKir dependence. It has been shown in the 
decays B° -> J/ip(KTr)f and B+ -> ir+(Kir)* ° Q that 
the amplitude behavior is consistent with that observed 
by LASS except for a constant phase shift. We allow 
an unconstrained overall shift, again with the index J 
suppressed, (So + ASq x Q) between the LASS amplitude 
phase and either the vector (J = 1) or the tensor (J = 2) 
resonance amplitude phase. 

The parameters £ describe the background or the re- 
maining signal PDFs. They are left free to vary in the fit 
for the combinatorial background or are fixed to the val- 
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FIG. 2: Projections onto the variables m& (a), m K ^ (b), 
AE (c), and tues (d) for the signal B° — > 4>(Kir) candi- 
dates. Data distributions are shown with a requirement on 
the signal-to-background probability ratio calculated with the 
plotted variable excluded. The solid (dashed) lines show the 
signal-plus-background (background) PDF projections. 



ues extracted from Monte Carlo (MC) simulation [12j and 
calibration B-decay channels for the exclusive B decays. 
We use a sum of Gaussian functions for the parameteri- 
zation of the signal PDFs for AE, rngg, and T . For the 
combinatorial background, we use polynomials, except 
for Toes an( i distributions which are parameterized by 
an empirical phase-space function and by Gaussian func- 
tions, respectively. Resonance production occurs in the 
background and is taken into account in the PDF. 

We observe a nonzero yield with more than 9a sig- 
nificance, including systematic uncertainties, in each of 
the three B° — > cf>K* decay modes. The significance 
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FIG. 3: Projections onto the variables Hi (a), H2 (b), Q $ 
(c), and the differences between the Q«I> projections for events 
with Hi Ti.2 > and with Hi H2 < (d) for the signal 
B° — ► 0A"*(892)° candidates following the solid (dashed) line 
dehnitions in Fig. [2] The D^-meson veto causes the sharp 
acceptance dips seen in (a). 
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TABLE I: Fit results for each mj& range and signal component: the reconstruction efficiency e r oco obtained from MC simulation; 
the total efficiency e, including the daughter branching fractions |3|; the number of signal events n s i g ; statistical significance 
(S) of the signal; the branching fraction B; and the flavor asymmetry Acp- The branching fraction B(B° —> ((i(K-k)q ) refers to 
the coherent sum |A res ± A non _ res | 2 of resonant and nonresonant J p — + Kir components pj and is quoted for rriKit < 1-6 GeV, 
while the B(B° —> 4>Kq (1430)°) is derived from it by integrating separately the Breit-Wigner formula of the resonant |v4 re s| 2 
Kir component [3] without rriKn restriction. The systematic errors are quoted last. 

Mode ma (GeV) model e rcco (%) e (%) n s j g (events) S (a) B (1CT 6 ) Acp 

4>K* (892)° 0.75 - 1.05 Breit-Wigner [3] 35.0 ±1.7 11.5 ± 0.6 406 ± 29 ± 15 21.0 9.2 ± 0.7 ± 0.6 -0.03 ± 0.07 ± 0.03 
(f)KZ(U30)° 1.13-1.53 Breit-Wigner [3] 27.1 ± 1.3 4.4 ± 0.2 133 ± 19 ± 7 9.7 7.8 ±1.1 ±0.6 -0.12 ± 0.14 ± 0.04 
(t>(Kir)* ° 1.13-1.53 LASS [7] 23.4 ±1.1 7.7 ± 0.3 147 ± 23 ± 7 9.8 5.0 ±0.8 ±0.3 ±0.17 ±0.15 ±0.03 

^o(1430)° Breit-Wigner [7] 4.6 ± 0.7 ± 0.6 



is defined as the square root of the change in 2 In £ 
when the yield is constrained to zero in the likelihood 
C. In Figs. [5HH we show projections onto the variables. 
In Tables HI and HT1 the rij, Aj, and £ = {/V, f±, <^>||, 
4>i_, So, A° CPl Aq P , A</>||, Acj)±, AS } parameters of the 
B° -> 0AT*(892) decay or the 0AT 2 *(143O)° and (j>(Kn)tf 
decays are obtained from the fit in the lower or higher 
ttik-k range, respectively. 

The nonresonant K + K~ contribution under the (f> is 
accounted for with the B° — ► frjK* category. Its yield 
is consistent with zero in the higher rriKix range and is 
89 ± 18 events in the lower ttik-e range. The uncertainties 
due to interference are estimated with the samples 

generated according to the observed K + K~ intensity and 
with various interference phases analogous to 5q in Kit. 
These are the dominant systematic errors for the £ pa- 
rameters of the B° -> 4>K*(892f decay. 

We vary those parameters in £ not used to model com- 
binatoric background within their uncertainties and de- 
rive the associated systematic errors. We allow for the 
flavor-dependent acceptance function and the reconstruc- 
tion efficiency in the study of asymmetries. The biases 
from the finite resolution of the angle measurement, the 
dilution due to the presence of fake combinations, or 
other imperfections in the signal PDF model are esti- 
mated with MC simulation. Additional systematic un- 
certainty originates from B background, where we esti- 





FIG. 4: Same as Fig.^a) and^b), but for the signal B° 
4>K% (1430)° and (f>(Kir)a° candidates combined. 



mate that only a few events can fake the signal. The 
systematic errors in efficiencies are dominated by those 
in particle identification and track finding. Other sys- 
tematic effects arise from event-selection criteria, </> and 
K*° branching fractions, and number of B mesons. 

In the lower ttikk range the yield of the 4>(Kit)q con- 
tribution is 601^ events with the statistical significance 
of 7.9cr, including the interference term. The dependence 
of the interference on the Kir invariant mass @, @| allows 
us to reject the other solution near (27r — (f)n,Tr — <fij_) 
relative to that in Table HJ for the B° -> ^*(892)° de- 
cay with significance of 5.4c, including systematic uncer- 
tainties. We also resolve this ambiguity with statistical 
significance of more than 4<r with the B° or B° decays 
independently. Because of the low significance of our 
measured f\\ = (1 — /l — fx) (2.9c) and fx (1.6a) in the 
B° — > 0i^|(143O)° decay we have insufficient informa- 
tion to constrain <f>n and <px at higher significance and to 
measure five asymmetries, and so we fix these asymmetry 
parameters to zero in the fit in the higher rriK n range. 

The (V — A) structure of the weak interactions and the 
s-quark spin flip suppression in the diagram in Fig. [T] (a) 
suggest \A \ > > [3,01. This expectation 

is consistent with our measurements in the vector-tensor 
B° — » ^i^l (1430)° decay, but disagrees with our observed 
vector- vector polarization. In the B° — > cf>K* (892)° de- 
cay we obtain the solution 0m ~ <j>x without discrete 
ambiguities. Combined with the approximate solution 
/l ~ 1/2 and fx — (1 — /l)/2, this results in the approx- 
imate decay amplitude hierarchy |Ao| — \A + \\ » 
(and |I | ~ » 

We find more than 5<r (4cr) deviation, including sys- 
tematic uncertainties, of <^±(0ii) from either ir or zero in 
the B° 4>K* (892)° decay, indicating the presence of 
final-state interactions (FSI) not accounted for in naive 
factorization. The effect of FSI is evident in the phase 
shift of the cosine distribution in Fig. [3] (d). 

Our measurements of eight CP-violating parameters 
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TABLE II: Summary of polarization results. The dominant 
fit correlation coefficients (C) are presented for the <f)K* (892)° 
mode where we show correlations of So with <f>u/tf>j. and of A<5o 
with A0||/A0 ± . For the 07^1(1430)° mode, the dominant 
values of C are 32% for (So,(f>\\) and 26% for (4>^,(f)±). 



h 
f± 



8 
A 

•Acp 

Acp 
A0|| 
A^x 
A5 



B° -> 0J<"2(143O) B°->^*(892) 

0.853+° ggg ± 0.036 0.506 ± 0.040 ± O.Of 5 

0.045to.o4o ± °- 013 0.227 ± 0.038 ± 0.013 
2.90 ± 0.39 ± 0.06 2.31 ± 0.14 ± 0.08 



H0.55 
-0.87 

,+0.12 



2.24 ±0.15 ±0.09 
2.78 ±0.17 ±0.09 
-0.03 ±0.08 ±0.02 
-0.03 ±0.16 ±0.05 
±0.24 ±0.14 ±0.08 
±0.19 ±0.15 ±0.08 
±0.21 ±0.17 ±0.08 



}-53% 

} 61% 
37%/27% 
}-51% 

} 61% 

37%/27% 



rule out a significant part of the physical region and are 
consistent with no CP-violation in this decay. Signifi- 
cant nonzero CP- violating parameters would indicate the 
presence of new amplitudes with different weak phases. 
The A(/>j_ and A0|| are particularly interesting due to 
sensitivity to the weak phases of the amplitudes with- 
out hadronic uncertainties 0, such as the relative weak 
phases of A + i and Aq, while the CP- violating A<5o pa- 
rameter represents potential differences of weak phases 
among decay modes. 

In summary, we have performed an amplitude analysis 
and searched for CP-violation in the angular distribu- 
tion with the B° — > <fiK*° decays with the tensor, vec- 
tor, and scalar K*°. Our results are summarized in Ta- 
bles U and [ill and supersede corresponding measurements 
in Ref. [J]. Our vector-tensor results are in agreement 
with quark spin flip suppression and Aq amplitude domi- 
nance, whereas the vector- vector mode contains substan- 
tial A + i amplitude from a presently unknown source ei- 
ther within or beyond the standard model Q. 
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